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ABSTRACT 

We present catalogues of faint 1 .4-GHz radio sources from extremely-deep Very Large Array 
pointings in the Lockman Hole, the Hubble Deep Field North (HDF-N) and ELAIS N2. Our 
analysis of the HDF-N data has produced maps that are significantly deeper than those pre- 
viously published and we have used these to search for counterparts to submm sources. For 
each of the fields we have derived normalised differential source counts and in the case of 
the HDF-N find no evidence for the previously-reported under-density of sources; our counts 
are entirely consistent with those found for the majority of other fields. The catalogues are 
available as an online supplement to this paper and the maps are also available for download. 

Key words: galaxies: starburst - galaxies; formation - cosmology; observations - cosmol- 
ogy; early Universe 



; 1 INTRODUCTION 

, Our knowledge of the evolution of galaxies and active galactic nu- 
clei (AGN) arguably owes as much to radio surveys as to those in 
any other band. Our knowledge of the cosmological evolution of 
radio-loud AGN is an obvious example (Dunlop & Peacock 1990), 
and one which foretold with remarkable accuracy the history of ac- 
cretion luminosity, and even that of star formation (Dunlop 1998). 

Recent work has shown that deep (/xJy) radio surveys have 
the power to probe star-forming galaxies, unbiased by dust obscu- 
_ ration, out to very high redshifts. Moreover, deep radio data have 
. proved to be extraordinarily complementary to datasets obtained 
' in the submillimetre (submm), infrared (IR) and X-ray wavebands 
' (Ivison et al. 2002; Donley et al. 2005). 

Unfortunately, the radio waveband has not benefited from the 
technological improvements and investment seen at shorter wave- 
lengths. This is in stark contrast with our mapping capability in the 
submm and IR wavebands, which has improved by many orders of 
magnitude in barely a decade. 

At 20 cm, the National Radio Astronomy Observatory's 
(NRAO's) Very Large Array (VLA) has been the dominant survey 
facility for several decades, with the Australia Telescope (ATCA), 
MERLIN and Westerbork (WSRT) synthesis arrays all having 
found significant niches. The most recent major radio surveys were 
the NRAO VLA Sky Survey (NVSS), covering 82 per cent of the 
sky to CT = 450 ^iJy beam"^ (Condon et al. 1998) and FIRST (Faint 
Images of the Radio Sky at Twenty cm, the radio equivalent of 
the Palomar Observatory Sky Survey) which mapped 10* deg'^ to 
around cr = 140 /ijy beam"^ (White et al. 1997). 

Going deeper, radio surveys in a plethora of fields have 
reached noise levels of a ~ 10-20 /ijy beam~^, the confusion limit 
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at 20 cm for arrays smaller than ~5 km. These surveys, often un- 
published but totalling around 10 deg^ , have been carried out with a 
variety of instruments (predominantly the VLA, but also the WSRT 
and ATCA) and include the Hubble Deep Field North (Richards 
2000), the Hubble Deep Field South (Huynh et al. 2005), the 
Spitzer First Look Survey (Condon et al. 2003; Morganti et al. 
2004), the VLA-VIRMOS Deep Field (Bondi et al. 2003) and the 
Phoenix Deep Survey (Hopkins et al. 2003). 

Given the rarity with which surveys have probed signifi- 
cantly below the a ~ lO^Jybeam^^ level (Ivison et al. 2002; 
Muxlow et al. 2005) it is natural to assume that to do so must 
require enormous allocations of observing time. This is not the 
case. Here, we describe three radio pointings which, together, re- 
quired less than one week of VLA time. Each covers an area of 
~320arcmin^ with 1.5-arcsec resolution, two of them to a noise 
level of ~5 /iJy beam~^, the other to ~10 ^Jy beam^^. 

Two of the VLA pointings lie within two regions of the 
sky which have very low Galactic backgrounds (<0.5 MJy sr^^ at 
100 ^m and E{B — V) < 0.01) and which were selected for the 
Spitzer Wide-area InfraRed Extragalactic (SWIRE) survey. These 
are the 'Lockman Hole' and ELAIS N2. The Lockman Hole is 
a 4 X 3-deg2 region centred at 10'M8"\ +57°04' (J2000) with 
~ 5 X 10^^ cm~^ of Hi (Lockman, Jahoda & McCammon 1986). 
This makes it uniquely well suited to panchromatic, deep studies 
of the Universe. The final VLA pointing is located in the Hubble 
Deep Field North (HDF-N), part of the Great Observatories Ori- 
gins Deep Surveys (GOODS, Dickinson et al. 2003). The ELAIS 
N2 and HDF-N data have been previously discussed by Ivison et al. 
(2002) and Richards (2000), respectively. A new reduction of the 
HDF-N 1.4-GHz VLA data (Morrison et al., in preparation) has 
also recently featured in Pope et al. (2006). 

We describe our observations in §2, the data reduction in §3, 
the extraction of sources in §4, present differential source counts 
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for all three fields in §5 and discuss our findings in §6. Throughout 
we adopt a cosmology, with fim = 0.3, = 0.7 and Ho = 
70kms"^Mpc"\ 



2 OBSERVATIONS 

Deep, high-resolution, wide-field radio images of three northern 
legacy fields were obtained using the NRAO's^ VLA. These were 
the '8-mJy Survey field' located within the Lockman Hole (Ivi- 
son et al. 2002), ELAIS N2 (Ivison et al. 2002) and the HDF-N 
(Richards 2000). 

For the 8-mJy Survey region, we used the VLA to acquire 
75 hr of data with a 4:1 ratio of time spent in the 'A' and 'B' con- 
figurations. For ELAIS N2, we acquired 20 hr of data, with a 3:1 
ratio of A:B. The relatively low integration time was due to an un- 
resolved ~120-mJy radio source, immediately evident in the field, 
and sufficiently bright that dynamic-range limitations (<5,000 for 
the VLA) would never justify a longer integration. For the HDF-N, 
we gathered all the available data from the archive, a total of 45 hr, 
all in 'A' configuration^ . Observing parameters for each field are 
summarised in Table 1 . 

Data were recorded every 5 s (3.3 s for the HDF-N) across 
two bands, these having lower-edge frequencies of 1365 and 
1435 MHz. The correlator was configured (mode '4') to produce 
seven 3.125-MHz channels per band in both right and left-circular 
polarisations (RCP & LCP). This is a commonly-used observing 
mode at low frequencies, the chosen bands being relatively free 
of radio-frequency interference (rfl) and the use of multiple chan- 
nels per band greatly reducing the effects of bandwidth smearing. 
During each day's observations, approximately 300 s were spent 
on 3C 48 (0137-1-331) for ELAIS N2, or 3C 286 (1331-1-305) for the 
Lockman Hole and the HDF-N; these data were used to set the flux 
density scale. The nearby, unresolved source, 1035-1-564, was used 
as the phase, amplitude and bandpass calibrator in the Lockman 
Hole, and was observed for 100 s every hour. Similarly, scans of 
1625-1-415 were used to calibrate the ELAIS N2 data, and 1313-1-675 
for the HDF-N. The remainder of the time was spent integrating on 
the target fields themselves, typically for 45-50 min between cali- 
bration scans. When 1331-1-305 was not observed, the flux density 
scale was set using boot-strapped flux densities for 1035-1-564. 



3 DATA REDUCTION 

The reduction of the data described in §2 followed a simple but 
time-intensive formula; all data were reduced using the NRAO 
MPS package. First, raw data were loaded using FILLM with 
DOUVCOMP=— 1 to minimise dynamic-range issues, with data 
weights calculated based on the nominal sensitivity and integration 
time. Corrupted data were identified visually and flagged using a 
variety of tasks. Flux calibration proceeded in the standard fashion 
using CALIB, antenna gain correction factors for all sources being 
derived from the flux calibrator data, 3C 48 or 3C 286, where ap- 
propriate. The bandpass shape was calibrated using the unresolved 



^ NRAO is operated by Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 

^ An additional 8 hr of data was contained in the archive, but as this was 
recorded several weeks before the bulk of the data and with a longer (10 s) 
integration time we did not include it in our analysis. 
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Figure 1. Wide-field contour map of the ELAIS field demonstrating the 
location of the central 37 facets that are used to map the VLA primary 
beam as well as the many other (smaller) facets that were simultaneously 
mapped in order to remove the response of sources beyond the primary 
beam. Contours outwith a facet are imaging artefacts and not real sources. 



calibrator, once per scan. The solutions so found were then inter- 
polated onto the target data, thus taking into account any time vari- 
ability of the antenna bandpasses. Finally, anteima-based amplitude 
and phase corrections were found for each calibrator scan and again 
interpolated onto the target data. 

Due to the very large size of the datasets and the associated 
computing times involved in their mapping and self-calibration, it 
was desirable to reduce the amount of data as much as possible 
without degrading the final maps. This was achieved using UBAVG 
to average data from different days together, the amount of time 
averaging set, per baseline, by the requirement that the peak flux 
density of sources at the edge of the field were reduced by less than 
1 per cent due to time-averaged smearing. For those datasets where 
more than one VLA configuration was present, each configuration 
was first averaged in this manner before being combined together 
into a single file. 

The calibrated data were imaged and deconvolved with 
IMAGR. In all cases a uniform data weighting scheme was used, 
with the ROBUST parameter set to zero. Making a single map of 
each, very large, field was impractical and so instead we followed 
the methodology of Owen et al. (2005), the primary beam (having 
an extent of ~0.5° for a 25-m antenna at 1.4 GHz) being approxi- 
mated with a total of 37 reasonably-sized (1024 X 1024 pixel) maps 
(or 'facets'). Additional maps were also made of sources detected 
beyond the main lobe of the primary beam, this being necessary 
in order to remove their sidelobes which would otherwise protrude 
into the area of the sky covered by the central facets; the clean 
components from these sources could then also be included in sub- 
sequent self-calibration. An illustration of this technique is given 
in Fig. 1 for the ELAIS field. Areas of CLEANing were restricted 
using boxes placed around sources, this also making it easier to use 
the CLEAN components when self-calibrating. This was performed 
using CALIB after each mapping run, solving ultimately for both 
amplitude and phase, but initially in phase only. The shortest solu- 
tion intervals (after averaging of data in ubavg) could be as low 
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Table 1. Observing dates, included VLA configurations, total integration times (target and calibrators), la rms noise and synthesised beam parameters (fwhm) 
for all fields. 



Field 


Observing dates 


Project code 


Configurations 


Integration time 


1(T rms noise 


Synthesised beam 










(hr) 


(juJybeam"'^) 


(arcsec^) 


8-mJy Survey Field 


2001 Jan - 2002 Mar 


AD432, AI088, AI098 


'A'.'B' 


75 


4.6 


1.26 X 1.32 


ELAIS N2 


2001 Jan - May 


AD432, AI091 


'A'.'B' 


20 


9.6 


1.48 X 1.45 


HDF-N 


1996 Nov -1996 Dec 


AR368 


'A 


45 


5.8 


1.52 X 1.51 



as ~1 min, but were generally longer for the 8-mJy Survey region 
(Lockman Hole) as there are fewer bright sources in this field. 

Mapping and self-calibration were performed separately on 
the RCP and LCP data for all fields (excepting ELAIS N2 which 
is the least sensitive of our datasets) in order to remove problems 
associated with the VLA 'beam squint' . At the point at which it was 
determined that further processing was not necessary, the RCP and 
LCP maps were combined to produce Stokes /. The final maps are 
shown in Fig. 2, where all 37 Stokes 7 facets have been combined 
into a single image of each field; in each case a model of the VLA 
primary beam has been used to correct the brightness of each pixel. 
Also shown are zooms of a 400 x 400-arcsec^ region of each map, 
with this area shown on the full-field maps as a box. 

We measured the Icr rms noise in the maps at various points 
in the maps away from bright sources and obvious image artefacts 
caused by residual calibration and deconvolution errors. Such prob- 
lems are well illustrated by the full primary-beam map of the HDF- 
N in Fig. 2. A number of sources can be seen that are clearly the 
source of enhanced noise, particularly the two extended sources in 
the southern half of the map. That near 12'' 37™, +61° 57' has a 
particularly complicated structure and the lack of short spacings in 
this 'A' configuration-only data almost certainly results in larger- 
scale extended structure being missed, making accurate decon- 
volution impossible. The second source, near 12'^ 37"", +61° 57' 
is the brightest source in the field, with a peak brightness of 
25.3 mJy beam"'^, as well as having one of the most complex struc- 
tures in the field. Another extended source located beyond the pri- 
mary beam towards the North-East is responsible for the resid- 
ual sidelobe structure in this quadrant of the map. However, for 
all fields the measured noise was close to the theoretical value 
as measured using the VLA exposure calculator^ . The rms noise 
for each field was: 8-mJy Survey (4.6 /iJy beam" '^), ELAIS N2 
(9.6 /iJy beam" ^) and HDF-N (5.8 ^Jy beam" ^). These measure- 
ments are also shown in Table 1 along with the synthesised beam 
parameters. 

A surprising feature of the maps for all three fields was 
the presence of sidelobes associated with the brighter sources 
(>1 mjybeam"') which proved to be undeconvolvable. These 
were also noted, for the HDF-N, by Richards (2000) who sug- 
gested that the cause was a correlator problem associated with 
the 3.3-s integration time. In common with these authors we find 
that the sidelobes are characterised by being radially orientated to- 
wards the phase centre of the observations and have an amplitude 
of ~ten per cent. In addition, the sidelobe pattern is noticeably 
stronger on the side of the source nearest the phase centre. An ex- 
ample of these sidelobes is shown in Fig. 3 for the brightest source 
in the HDF-N. The sidelobes only appear to affect the maps in an 
area immediately around each source, the theoretical rms noise be- 



^ http://www.vla.nrao.edu/astro/guides/exposure/ 
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Figure 3. Undeconvolvable sidelobes around the brightest source in the 
HDF-N; these point directly towards the phase centre of the observations. 
The greyscale is given in injybeam"^ and the synthesised beam (fwhm) 
is illustrated at bottom-left. Axes give the offset in arcsec from the centre of 
the facet. 



ing achieved elsewhere in the maps, and therefore we conclude that 
their effects on the results are not particularly pernicious. 



4 SOURCE EXTRACTION 

The process of extracting the position, flux density and shape/size 
of each source in the data was designed to be as automated as 
possible. An initial source list was first formed using the AiPS 
task SAD. This searches for areas in maps which lie above a cer- 
tain threshold ('islands') and fits one or more (up to a maximum 
of four) two-dimensional Gaussians to these. It subsequently esti- 
mates the source size and orientation by deconvolving the CLEAN 
beam from the Gaussian fit. sad was not run on the maps produced 
by IMAGR but on signal-to-noise ratio (SNR) maps which were cre- 
ated by dividing each IMAGR map with its noise map, this having 
been constructed using the AIPS task RMSD. In this manner it was 
possible to restrict SAD to search for only those sources that had a 
peak flux density in excess of 4(t, where a is the local rms noise. 
The noise in the maps, whilst for the most part uniform, does vary 
with position, mainly due to residual calibration and deconvolution 
errors as well as the afore-mentioned radial sidelobe pattern. 
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Figure 2. Stokes I maps of of the 8-mJy Survey field (top), ELAIS N2 field (middle) and the Hubble Deep Field North (bottom). Right: full mosaic, with a 
diameter of 45 arcmin. Left: small 400 x 400-arcsec portion of the area boxed in the panels on the right. Grey scales are plotted from 0-100 juJy beam~i . 
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Table 2. Extract (first ten sources) from the HDF-N catalogue, the full listing of which (along with the 8-mJy and ELAIS N2 fields) is available online. 
Quantities in brackets give the Icr error. The third column gives the distance of the source from the phase centre. Sources are listed in order of Right 
Ascension. Unresolved sources have no entries in the size columns and have integrated flux densities equal to the peak flux density. The final column includes 
a note signifying whether the source is so extended that its flux densities have been estimated using tvstat ('E'), or whether two Gaussians have been grouped 
together to form a single source ('D'). 



Right Ascension 


Declination 

(-"") 


Radial distance 
(arcmiii) 


Peak flux 
(//.lybcani 


Total flux 
(//.ly) 


Deconvolved source 
Major (mas) Minor (mas 


size Note 
1 P.,\. ( -') 


12:33:50.552 (0.027) 


62:19:35.36 (0.17) 


21.8 


333 (52) 


434 (89) 


1.9 


0.0 


86.1 


12:33:51.801 (0.004) 


62:16:07.57 (0.01) 


20.9 


2206 (45) 










12:34:00.031 (0.035) 


62:12:22.08 (0.12) 


19.7 


192 (36) 










12:34:03.510(0.009) 


62:14:20.52 (0.03) 


19.4 


637 (33) 










12:34:09.917 (0.033) 


62:03:57.34 (0.15) 


20.7 


237 (45) 










12:34:10.672(0.015) 


62:06:15.76(0.06) 


19.7 


441 (36) 


461 (53) 


0.7 


0.3 


81.1 


12:34:11.238 


62:16:16.78 


18.7 


2349 


7642 


10.4 




126.5 E 


12:34:14.475 (0.024) 


62:21:09.29 (0.13) 


19.8 


238 (35) 


240 (50) 


0.9 


0.0 


85.1 


12:34:16.104(0.018) 


62:07:22.18 (0.08) 


18.7 


307 (30) 


365 (49) 


1.6 


0.0 


46.0 


12:34:16.532(0.032) 


62:21:41.77(0.15) 


19.8 


195 (36) 











A small 'postage-stamp' map was produced at the position of 
each fitted model component, with the model Gaussian itself also 
plotted. By visually inspecting each map it was possible to remove 
components that obviously did not correspond to a real source, but 
instead to a deconvolution artefact or sidelobe. In this first stage of 
model fitting, many components were also fitted to large extended 
sources whose size, shape, flux and position could not be reason- 
ably approximated by one or more Gaussians. With both of these 
component types removed from further consideration in the auto- 
mated source extraction process, the remaining components were 
model fitted again, this time using jmfit. 

From this point, close pairs (where the sum of the major axes 
was greater than or equal to the inter-component distance) were 
fitted simultaneously. The area made available to jmht was cho- 
sen per sotirce, again based on the component sizes as well as the 
inter-component distances where appropriate. Using the geomet- 
ric mean of the major and minor axes, those components where 
the undeconvolved component size was smaller than the restoring 
beam (corrected for bandwidth smearing) were considered to be 
uru-esolved and their sizes held fixed at the size of the (smeared) 
restoring beam. This is the case for approximately half the sources 
that appear in the final catalogues, for each field. 

The results of this model fitting were then used to reject all 
source components with a peak flux density below 5a. A final 
round of model fitting was then performed, this time on the total 
intensity maps (not the SNR maps) with corrections made both for 
the attenuation of the primary beam and the effect of bandwidth 
smearing. At all times the model fitting was performed on the indi- 
vidual map facets and not on the combination maps shown in Fig. 2. 
Due to this there are occasions when a source close to the map edge 
is detected in more than one facet. If so, the fit from the map where 
the source Ues closest to the centre of the map is retained and the 
others discarded. 

For those sources that could not be represented by Gaussians 
(18, 10 and 7 in the ELAIS, Lockman Hole and HDFN maps re- 
spectively), their total flux densities were measured from the maps 
by stmmiing all the pixels adjudged, by eye, to make up the source 
(using TVSTAT). The position of such a sotirce (marked with an 'E' 
in the catalogue) corresponds to the peak, the catalogues also giving 
the largest angular size of the source and its orientation (measured 
using IMDIST). Of the sources fitted with Gaussians, those com- 
ponents which have similar integrated flux density ratios (within 



a factor of two) and which are so closely separated that the com- 
ponents are clearly merging are considered to be the same source. 
There were also occasions when more widely separated compo- 
nents were clearly associated due to a combination of their rela- 
tively small proximity given their high brightness and/or the pres- 
ence of faint extended emission between and on a Une joining the 
two components. In all cases these were double sources and as with 
the extended sources their catalogue entries give the position of 
the brighter component and arc marked with a 'D'. The number 
of sources which have been formed from multiple components is 
small: four in the HDFN, four in ELAIS and none in the Lockman 
Hole. 

Finally, a number of checks for missing sources were made, 
both by inspecting the residual maps made by SAD and by making 
maps of each facet with the positions of the sources plotted. 

The complete catalogues are available as on online supple- 
ment to this paper and we illustrate the catalogue format for the first 
ten sources of the HDF-N in Table 2. The total number of sources in 
each catalogue is 537 (HDFN), 200 (ELAIS) and 506 (Lockman). 
We are also making our maps available to the community and these 
can be downloaded, in FITS format, from the internet*. 



5 SOURCE COUNTS 

We have used our source catalogues to derive the differential nor- 
malised source counts for the three fields. In doing this it is cru- 
cial to correct for a number of effects that lead to a reduction in 
the number of sources detected, particularly at low flux densities. 
These include: 

(i) Increase in map noise due to attenuation of primary beam. 

(ii) Reduction in source peak flux density due to bandwidth 
smearing. 

(iii) Resolution bias. Sources with integrated flux densities > 
5(7 are missed if they are sufficiently extended that their peak flux 
density < Stx. 

(iv) Reduced reliabiUty of SAD at recovering sources when the 
signal to noise is low. 

ftp://ftp.roe.ac.uk/pub/adb/vla.deep.fields/ 
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Figure 4. Cumulative distribution of catalogued ratios of peak and inte- 
grated flux densities (Sp/Si). For each field it can be seen that approxi- 
mately half of the sources have identically equal integrated and peak flux 
densities (fitted as unresolved sources in jmht.) 



In order to estimate the magnitude of these biases we have 
simulated a radio population with the same characteristics as those 
that form our catalogues. We insert 50 sources into each of the 37 
residual maps output by sad, after these have been corrected for 
the attenuation of the primary beam and the effect of bandwidth 
smearing. The source integrated flux densities, Si, are chosen ran- 
domly from a logaritttmic distribution; the lowest flux density is set 
a little below the ~5cr flux Umit of the field (20 /uJy beam~^ for 
the HDFN and Lockman Hole and 40 pJy beam~^ for ELAIS) and 
the highest at a level sufficient to encompass the brightest sources 
used in our source counts. 

We model our source sizes from the data themselves. The area 
of each modelled source is related to the ratio of the peak flux den- 
sity, Sp, and Si. Having formed the cumulative distribution of peak 
over integrated flux density for each field (Fig. 4), we pick a value 
at random and use this to calculate a modelled source size, Omod, 
using 



Si ^ 
Sp 



1 1/2 



(1) 



where Omaj and ( 



Jmm are the major and minor axes of the restoring 
beam. All modelled sources are circular. 

Once the flux density and size parameters are determined, 
sources can be added to the residual maps (after convolution with 
the restoring beam) using IMMOD, the positions of each source also 
being chosen at random. Our source extraction sequence is then run 
on the simulated maps, exactly as with the observed data, and a note 
made as to whether a source was recovered or not. 

The source counts themselves are formed by counting the 
number of observed sources, A'^, within bins of integrated flux den- 
sity which start at 5a and increase by factors of 1.5. For each bin, 
the smallest distance from the phase centre, R, necessary for all 
sources within that bin to be included is noted* . The somce density 
(Jy~^ sr~^) is then calculated for each bin using the area corre- 



Due to increasingly large errors in the polynomial that describes the pri- 
mary beam correction, we do not include any sources beyond 20 arcmin in 
our source count calculations. 



spending to the radius R and then Euclidean-normalised by multi- 
plying by the average flux density of all the sources in the bin, S, 
to the power 2.5. The error on each bin's value is assumed to be 
Poissonian i.e. VN. Only bins that contain at least nine sources are 
included. 

Each bin's count is then corrected for the biases described 
above: for all those simulated sources which fall into the bin, the 
fraction that are recovered (over the same area corresponding to the 
bin's radius R) is calculated and its inverse used as a multiplicative 
corrective factor, C. The bin error is also scaled by this same factor. 
The final source counts are given in Tables 3, 4 and 5 and plotted 
in Fig. 5 (including the HDF-N source counts of Richards.). 



6 DISCUSSION 

Fig. 5 only significantly samples the so-called 'upturn' in the nor- 
malised differential source counts, the point at flux densities below 
approximately 1 mJy where the steep decline seen at higher flux 
densities, just visible in our plot, is arrested. We find that there is 
general agreement between the source counts for each of the three 
fields and that these are consistent with the results found by other 
authors for other fields; the HDF-N source counts, for example, are 
extremely similar to those found by Seymour et al. (2004) for the 
O*" XMM-Newton/ROSAT fisld. The counts for the Lockman Hole 
are lower, especially for Si < 200/Jy, but this is not altogether 
surprising given that less sources were detected here than in the 
HDF-N, despite the Lockman Hole data being more sensitive. 

The correction factors quoted in Table 3 are higher than those 
quoted by some other authors, but it is important to compare like 
with like. For instance, the correction factors of Richards (2000) 
do not include the attenuation of the primary beam or bandwidth 
smearing, this instead being factored into their quoted value of 
dN/ dS. Seymour et al. (2004), do include primary beam and chro- 
matic aberration effects in their correction factor and for their low- 
est bin this is equal to 2.56, similar to our corresponding value for 
the HDF-N. In general, source counts derived from single pointings 
of a radio telescope (e.g. Richards 2000; Seymour et al. 2004) will 
require large corrections to the lowest bins because the noise is so 
non-uniform over the beam area. Observations which mosaic sev- 
eral pointings together (e.g. Hopkins et al. 1998; Bondi et al. 2003) 
have more uniform noise and consequently lower corrections. 

Most noteworthy in Fig. 5 is the striking inconsistency be- 
tween the HDF-N source counts presented in this paper and by 
Richards (2000), the new counts being significantly higher. 



6.1 Disparity with the Richards (2000) results 

The finding that our source counts for the HDF-N are higher than 
those published by Richards (2000) is particularly interesting as 
the previous result has been noted for many years (beginning with 
Richards) as being anomalously low; the new counts on the other 
hand are extremely consistent with the results from other fields. 
This in itself would suggest that the source counts presented here 
are more accurate, but, given the disparity, it is important to perform 
a careful comparison with the original. 

In terms of number of sources, our catalogue contains many 
more, 537 versus 371. This is partly due to the new catalogue cov- 
ering a larger area; out to the same radius (20 arcmin), Richards 
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Table 3. Differential normalised source counts for the HDF-N. Tabulated are the bin upper and lower limits (Si and Su), average flux density (S), number of 
sources per bin {N), numbers of sources per bin from the Richards (2000) catalogue (ATj^qo). radius (in arcmin) out to which sources were counted (Ft), bin 
correction factor (C) and the normaUsed differential source covmts themselves. 



Si Su S N ATroo R C S'^-^dN/dS 

(/iJy) (/iJy) (/iJy) (arcmin) (Jy-*^'^ sr"-') 



30.0 


45.0 


40.2 


43 


1 


9 


2.50 


3.41 ± 0.52 


45.0 


67.5 


56.4 


105 


63 


13 


1.41 


3.50 ± 0.34 


67.5 


101.2 


80.9 


99 


96 


16 


1.56 


3.96 ± 0.40 


101.2 


151.9 


122.4 


94 


67 


19 


1.35 


4.34 ± 0.45 


151.9 


227.8 


184.5 


63 


59 


20 


1.14 


4.10 ± 0.52 


227.8 


341.7 


266.0 


42 


40 


20 


1.04 


4.17 ± 0.64 


341.7 


512.6 


418.8 


20 


22 


20 


1.00 


3.95 ± 0.88 


512.6 


768.9 


629.5 


11 


14 


20 


1.00 


4.01 ± 1.21 


768.9 


1153.3 


988.1 


9 


7 


20 


1.01 


6.83 ± 2.28 



Table 4. Differential normahsed source counts for the 8-mJy Survey Region of the Lockman Hole. 



Si Su S N R C S'^-^dN/dS 

(/iJy) (/iJy) (/iJy) (arcmin) (Jy^'^sr"^) 



23.0 


34.5 


29.3 


55 


9 


3.03 


3.14 ± 0.42 


34.5 


51.8 


42.0 


102 


11 


1.28 


2.70 ± 0.27 


51.8 


77.6 


63.4 


82 


15 


1.37 


2.32 ± 0.26 


n.6 


116.4 


98.0 


73 


18 


1.23 


2.56 ± 0.30 


116.4 


174.7 


143.4 


47 


20 


1.10 


2.05 ± 0.30 


174.7 


262.0 


220.0 


46 


20 


1.03 


3.67 ± 0.54 


262.0 


393.0 


309.4 


29 


20 


1.03 


3.62 ± 0.67 


393.0 


589.5 


490.6 


14 


20 


1.01 


3.61 ± 0.96 


589.5 


884.2 


696.6 


9 


20 


1.00 


3.68 ± 1.23 


884.2 


1326.3 


1099.5 


13 


20 


1.00 


11.09 ±3.07 



finds 370 sources'' compared to our 506. The remaining surplus is 
due to our maps being significantly deeper and shouldn't contribute 
to the discrepancy in the source coimts. 

We note that there are 34 sources in the Richards catalogue 
that do not appear in ours. We have scrutinised our maps at the po- 
sitions of these sources and in the majority of cases conclude that 
weak emission is present, but below our 5a peak cutoff. Most of 
these sources have a quoted SNR of between five and seven al- 
though one, with a total flux density of 1.9 mJy (SNR > 40) is 
certainly not real, there being blank sky at the quoted position both 
in our map and in FIRST. Different data weighting is one possi- 
ble explanation for why the other sources do not make it into our 
catalogue. 

By identically binning the sources from both catalogues we 
can directly compare the source coimts from each. We have binned 
the Richards sources in the same manner as our own and find, not 
surprisingly given the better sensitivity of the new maps, that in the 
lowest bins we detect many more sources (Table 3). In the subse- 
quent bins there is a general trend for our catalogue to contribute 
more sources, but in two cases the Richards catalogue has more. 
In general, the difference in the nimibers detected is not sufficient 
to accotmt for the systematic offset. For example, at ~180/iJy, 



The first source in the Richards catalogue is actually 28.3 arcmin from 
the phase centre. 



both catalogues have similar numbers of sources, but there is al- 
most a factor of two difference between the source counts. We have 
also calculated the integrated fiux density ratios for the 334 sources 
common to both catalogues and find (after removal of three outliers 
from Richards whose flux densities are in error by factors ^ 10) a 
mean of 1.02 and an rms of ~20%. 

Ultimately, we are unable to identify the cause of the lower 
HDF-N source counts. However, we note that it does not seem pos- 
sible to arrive at the values of n/uo from the information given 
in Table 3 of Richards. These should be obtainable by multiplying 
dN/dS by the correction factor and then by the average flux den- 
sity to the power 2.5. If we assume that column 2 (labelled (51.4)) 
is the appropriate measure of flux density, we find that the numbers 
obtained in this manner do not agree with and are in fact higher 
than those reported in Table 3 (apart from the highest flux density 
bin) and therefore in better agreement with our counts. 



6.2 New submm counterparts in the HDF-N 

An invaluable application of deep radio catalogues to cosmology 
has been their usefulness in identifying the optical counterparts to 
submm sources detected by e.g. the SubmilUmetre Common User 
Bolometer Array (SCUBA, Holland et al. 1999). The much smaller 
beam of the VLA at 1.4 GHz (1-2 arcsec) compared to that of 
SCUBA (~15 arcsec at 850 jjm) and the fact that submm sources 
often have associated radio emission makes it possible to get much 



(£) 0000 RAS, MNRAS 000, 000-000 



8 Biggs & Ivison 

Table 5. Differential normalised source counts for ELAIS N2. 



Si 


Su 


S 


AT 


R 


C 


52-5 dN/dS 


(MJy) 


(MJy) 


(A'Jy) 




(arcmin) 




(Jyi-5 sr-i) 


50.0 


75.0 


68.1 


13 


7 


3.03 


4.63 ± 1.28 


75.0 


112.5 


95.0 


29 


12 


1.67 


2.96 ± 0.55 


112.5 


168.8 


140.2 


30 


14 


1.30 


3.09 ± 0.56 


168.8 


253.1 


209.0 


29 


18 


1.41 


3.55 ± 0.66 


253.1 


379.7 


297.2 


19 


19 


1.16 


2.77 ± 0.64 


379.7 


569.5 


450.6 


14 


19 


1.06 


3.52 ±0.94 


569.5 


854.3 


702.9 


12 


20 


1.04 


5.41 ± 1.56 
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Figure 5. Normalised differential source counts for the HDF-N (open squares), the Lockman Hole (open diamonds) and ELAIS N2 (open circles). For 
comparison we have also plotted the Richards (2000) HDF-N source covmts (filled squares). 



more accurate positions for SCUBA sources which can then be 
used to identify the galaxy in the optical and infrared (e.g. Ivison 

et al. 2002). 

Interest in the HDF-N remains high, especially with the ad- 
vent of GOODS (Dickinson et al. 2003), a deep multi-wavelength 
survey which, amongst other things, seeks to constrain theories of 
galaxy evolution. Its northern field, GOODS-N, coincides with the 
HDF-N. Pope et al. (2005) present a recent SCUBA 850-^m cat- 
alogue of the GOODS-N field which contains 40 sources, several 
of which are previously unpublished. Of these, 27 have an optical 
coimterpart and twelve of these do not have known radio detections. 
In Fig. 6, a 20 X 20 arcsec^ VLA map centred on the submm posi- 
tion of each of these twelve sources is plotted. A circle indicates the 
2(7 positional error (SNR dependent, Ivison et al., in preparation) 
and a cross the position of the optical/near-IR identification. 



Two of the submm sources (GN03 and GN20) have a clear 
5a radio detection within the error circle and both of these are in- 
cluded in our HDF-N catalogue. Of particular interest is GN20, 
this previously being the brightest (21.3 mjy) submm source with- 
out an identified cotmterpart in a different waveband. It is, how- 
ever, clearly detected in our map, is slightly resolved and lies only 
1.3 arcsec away from the submm position. In the case of GN03, 
the radio/submm positional offset is 2.4 arcsec. Positions and flux 
densities for these two sources are given in Table 6. Similar flux 
densities for these sources are found by Pope et al. (2006). A third 
source, GN32, also has a radio source within the 2cr error circle, 
but this is already present in the Richards (2000) catalogue. 
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Figure 6. VLA 1.4-GHz maps centred on the positions of the 12 submm sources in the HDF-N described by Pope et al. (2005) as having an optical/near-IR, 
but no radio, counterpart. Each map is 20 x 20 arcsec-^ with contours plotted at —3, 3, 4, 5, 6 times the la rms noise. The circle describes the 2a submm 
position error based on the signal to noise ratio (Ivison et al., in preparation). The cross marks the position of the optical/near-IR coimterpart (Pope et al. 2005). 
EUipses in the bottom-left comers represent the telescope beam (fwhm). 
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Table 6. Submm sources with new radio counterparts identified in the new reduction of the HDF-N data. Submm names and positions are taken from Pope 
et al. (2005). 



Submm name Submm position Radio position Radio flux density (/iJy) 

GN03(SMMJ123608+621147) 12''36™08!9, +62° 12'53" 12'i36™08!665, +62°12'51'.'26 34±7 
GN20(SMMJ12371 1+622212) 12''37™11!7, +62°22'12" 12''37™11?875, +62°22'11'.'54 58+16 



7 CONCLUSIONS 

In this paper we have presented radio maps and catalogues from 
three legacy fields, the HDF-N, the 8-mJy Survey Region of the 
Lockman Hole and ELAIS N2. These are areas of the sky which 
continue to be observed at all wavelengths and with a variety of 
telescopes, resolutions, etc; we hope that by making our radio maps 
available, the continued exploitation of these fields will be made 
much easier given that the protracted tasks of assembly, calibration 
and mapping of the data will be not be necessary. 

The data presented here have been discussed before by other 
authors, but the new reduction of the HDF-N data is a marked 
improvement on that already published (Richards 2000) and has 
yielded a new source catalogue with more and fainter sources. This 
has been used to find radio counterparts to two submm sources 
which previously had none; both of these identifications have re- 
cently also been made by Pope et al. (2006) from a separate re- 
analysis of the HDF-N data (Morrison et al., in preparation). By 
going deeper still it may be possible to identify the radio counter- 
parts to those submm sources in the HDF-N which remain unde- 
tected, sources which otherwise might be considered spurious (e.g. 
Greve et al. 2004). Deeper, more reliable radio maps of the HDF-N 
should be available in the near future due to further VLA obser- 
vations at 1.4 GHz in the 'A, 'B', 'C' and 'D' configurations (P.I. 
Gleim Morrison). 

Source counts have been presented for each of our fields, down 
to a flux density limit of 5a. These are broadly consistent with 
one another, but the observed differences (such as the low counts 
seen in the Lockman Hole) might be a consequence of cosmic vari- 
ance. Our most important conclusion is that we find no evidence for 
the under-density of sources in the HDF-N reported by Richards 
(2000), instead obtaining results which are entirely consistent with 
those found for other fields. It has long been noted that the HDF- 
N source counts were anomalously low and this finding appears to 
have been spurious. The recent resurgence in interest in the HDF-N 
radio data, both in the re-analyses of the archived data by ourselves 
and Morrison et al. (in preparation) and the on-going effort to add 
new data, should confirm our results and allow the source coimts to 
be extended to still fainter flux densities. 



REFERENCES 

Bondi M. et al., 2003, A&A, 403, 857 

Condon J. J., Cotton W. D., Greisen E. W., Yin Q. F., Perley R. A., 
Taylor G. B., Broderick J. J., 1998, AJ, 1 15, 1693 

Condon J. J., Cotton W. D., Yin Q. F., Shupe D. L., Storrie- 
Lombardi L. J., Helou G., Soifer B. T., Werner M. W., 2003, 
AJ, 125,2411 

Dickinson M., Giavalisco M., The Goods Team , 2003, in The 

Mass of Galaxies at Low and High Redshift, p. 324 
Donley J. L. et al., 2005, AJ, 129, 220 

Dunlop J. S., 1998, in ASSL Vol. 226: Observational Cosmology 

with the New Radio Surveys, p. 157 
Dunlop J. S., Peacock J. A., 1990, MNRAS, 247, 19 
Greve T. R., Ivison R. J., Bertoldi F., Stevens J. A., Dunlop J. S., 

Lutz D., Carilli C. L., 2004, MNRAS, 354, 779 
Holland W. S. et al., 1999, MNRAS, 303, 659 
Hopkins A. M., Afonso J., Chan B., Cram L. E., Georgakakis A., 

Mobasher B., 2003, AJ, 125, 465 
Hopkins A. M., Mobasher B., Cram L., Rowan-Robinson M., 

1998, MNRAS, 296, 839 
Huynh M. T, Jackson C. A., Norris R. P., Prandoni 1., 2005, AJ, 

130, 1373 

Ivison R. J. et al., 2002, MNRAS, 337, 1 

Lockman F. J., Jahoda K., McCammon D., 1986, ApJ, 302, 432 

Morganti R., Garrett M. A., Chapman S., Baan W., Helou G., 

Soifer T., 2004, A&A, 424, 371 
Muxlow T. W. B. et al., 2005, MNRAS, 358, 1159 
Owen F. N., Keel W. C, Ledlow M. J., Morrison G. E., Windhorst 

R. A., 2005, AJ, 129, 26 
Pope A., Borys C, Scott D., Conselice C, Dickinson M., 

Mobasher B., 2005, MNRAS, 358, 149 
Pope A. et al., 2006, MNRAS, accepted (astro-ph/0605573) 
Richards E. A., 2000, ApJ, 533, 611 

Seymour N., McHardy 1. M., Gunn K. F, 2004, MNRAS, 352, 
131 

White R. L., Becker R. H., Helfand D. J., Gregg M. D., 1997, ApJ, 
475, 479 



ACKNOWLEDGMENTS 

ADB would like to thank Nuria Lorente for her help in getting var- 
ious perl scripts and IDL programs working and Nick Seymour for 
his thoughts on source counts. RJI acknowledges kind and patient 
help over the years from Ernie Seaquist, Frazer Owen, Bob Becker, 
Chris Carilli and Rick White. We thank the anonymous referee for 
many suggestions which greatly improved this manuscript. 



© 0000 RAS, MNRAS 000, 000-000 



